Fluctuating particle motion during shear induced granular compaction 
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Using a refractive index matching method, we investigate the trajectories of particles in a three 
dimensional granular packing submitted to cyclic shear deformation. The particle motion observed 
during compaction is not diffusive but exhibits transient cage effect, similar to the one observed in 
colloidal glasses. We precisely study the statistics of the step size between two successive cycles and 
observe that it is proportional to the shear amplitude. The link between the microscopic observations 
and the macroscopic evolution of the volume fraction during compaction is discussed. 

PACS numbers: 45.70.Cc, 05.40.Fb 



Granular material are often described as athermal sys- 
tems [liy. For grains greater than few microns, thermal 
fluctuations are negligible and play no role in the dy- 
namics. A granular assembly can then be trapped in a 
metastable configuration and will not evolve unless an ex- 
ternal perturbation such as vibration or shear is applied. 
The compaction of an initially loose granular sample rep- 
resents a simple configuration to study the dynamics of 
athermal system. 

Two configurations have been experimentally inves- 
tigated to study granular compaction. The first one 
consists in imposing vertical taps to a granular packing 
d • The second consists in imposing a periodic 
shear deformation to a sample contained in a box 0, • 
In both cases, when the amplitude of excitation is con- 
stant (the tap intensity or the shear amplitude), a slow 
logarithmic relaxation towards more compact states is 
observed 0, 0, 0- In contrast to this slow dynamics, 
rapid changes in volume fraction can be observed when a 
sudden change in excitation is imposed • The pack- 
ing dilates (resp. compacts) when the excitation ampli- 
tude increases (resp. decreases). Moreover the system 
exhibits a short term memory, i.e. two packings at the 
same volume fraction but being prepared differently do 
not evolve the same way when submitted to the same 
excitation. The coexistence of a slow and a rapid dy- 
namics and the observation of memory effects have moti- 
vated number of numerical or theoretical models mainly 
based on free volume concept and geometrical constraints 
H GH El El El G3 However, to our knowledge, no 
link has been made between the macroscopic dynamics 
and the microscopic structure of the packing during com- 
paction. Does the slow and rapid dynamics correspond 
to different grains rearrangements? What does the mem- 
ory effect mean for the particles position? In this study 
we address these questions by analyzing the individual 
motion of the grains during a cyclic shear compaction. 

In thermal system like colloidal suspensions, the micro- 
scopic motion of the particles has been recently studied in 
detail. Using confocal microscopy, Weeks et al. 0, 
have been able to follow the three dimensional trajecto- 



ries of the particles in colloidal hard spheres suspensions. 
They have shown that for dense systems near the glass 
transition the particles spend most of the time trapped 
in localized regions or "cages" and occasionally exhibit 
longer excursions. Considering the analogy often pro- 
posed between granular compaction and aging in glassy 
systems |2j , one can legitimately wonder if similar micro- 
scopic behaviors are present in granular compaction. 




FIG. 1: a) Experimental setup, b) Picture of the packing 
full of index matching liquid with 12 colored tracers. The 
fuzziness allows a better accuracy for particle tracking, c) 
Picture obtained with the laser sheet method. 

The experimental setup is described in details in 
and is presented in Fig. 1. Spherical glass beads (d = 3 
mm in diameter) contained in a parallelepipedic box (7.7 
cm wide, 10.2 cm deep and typically 10.5 cm high) are 
submitted to a periodic shear by inclining the two side 
walls. The top plate which applies a constant pressure 
on the packing is free to move vertically only. During 
compaction the top plate slowly goes down. The measure 
of its vertical position gives the bulk volume fraction. All 
the experiments are carried out in a quasistatic regime 
and the only control parameter is the amplitude of the 
cyclic shear measured by the maximum inclination of 
the side walls. 

The particles are tracked during compaction using an 
index matching method. The material is immersed in a 
mixture of 67% of turpentine and 33% of methylnaph- 
talene matching the glass refractive index of the beads. 
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FIG. 2: Compaction for = 5.4°. a) Volume fraction as 
function of the number of cycles, b) Examples of trajectories 
during 15000 steps. The disks give the beads size and indicate 
the initial position of the tracers, c) Examples of cages (tra- 
jectories ploted for time slots between 2500 and 5000 steps). 



represents the stroboscopic trajectory i.e. the successive 
positions of the particles measured after each shear cycle. 
Trajectories are first characterized by an overall decrease 
of the vertical coordinate as the number of cycles in- 
creases, reminiscent of the macroscopic compaction. The 
corresponding evolution of the volume fraction is plot- 
ted in Fig 2a. On top of this mean vertical displace- 
ment, one observes fluctuating motion characterized by 
ball-like regions as shown in the close-up Fig 2c. This 
could reveal caging process. The random motion of the 
particles seems to be trapped in a finite volume before 
escaping and being trapped again in another cage. No- 
tice that the change of cage corresponds to a displace- 
ment quite smaller than the particle diameter (about one 
fourth). When compaction goes on, the time spent in a 
cage becomes longer and longer and the particles even- 
tually move around a fixed position. The cage effect is 
reminiscent of what has been observed in colloidal glasses 
by Weeks et al. [01 • The main difference here is that 
the system is not time invariant and evolves towards more 
and more compact states. Cage properties changing with 
time and particle experiencing only few cage changes be- 
fore being trapped in its final location (between 2 and 5), 
we have not been able to get significant statistics about 
the cage sizes. 



The packing is then transparent and it is possible to fol- 
low few colored glass beads initially introduced in the 
packing (Fig. lb). The experimental procedure is the 
following. The top plate being removed, the particles 
are gently poured in the container filled with the index 
matching liquid. Colored particles are then carefully in- 
troduced in the packing and the top plate is slowly put 
into contact with the beads. The mean initial volume 
fraction is typically 0.61. Once the packing is ready, the 
shear deformation is imposed: the side walls are succes- 
sively inclined to an angle +6 then — 0, then put back at 
the vertical position. After each cycle, we measure the 
volume fraction from the position of the top plate and we 
take a front picture of the packing. The tracers locations 
in the x-y plan are determined from the picture using 
a particle tracking algorithm. The location accuracy is 
±15/ira. 

In order to observe trajectory correlations between ad- 
jacent particles, another visualization technics is devel- 
oped. Fluorescent dye is added to the mixture (with some 
alcohol to help dilution) and the packing is lighted with 
a vertical thin laser sheet in the x-y plane. One then ob- 
tains a picture where the particles appear as black disks 
(Fig lc). This technics allows to study short time corre- 
lation motion but does not allow to follow particles in the 
long time, as they can move perpendicular to the laser 
plane and eventually leave it. 

The first experiments are performed at a constant 
shear amplitude 0. An example of the particle motion 
obtained for = 5.4° is presented in Fig. 2b. The plot 
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FIG. 3: a) Normalized distribution P(Ar) computed at dif- 
ferent stages during compaction; a is the standard devia- 
tion. = 8.1°; (o) distribution obtained for 1 < n < 1000 
(a = 72/xra); (+) 4000 < n < 5000 (a = 54/xra); (x) 
9000 < n < 10000 (a = 52/xra). Solid line is the gaussian 
distribution, b) x as a function of Ar (see text). The line is 
a linear fit. 

However, in order to investigate how the random mo- 
tion evolves with compaction, we study the statistics of 
the displacement Ar = \J Ax 2 + Ay 2 experienced by the 
particles between two successive cycles. The distribution 
of P(Ar) obtained at various stages of the compaction 
process are plotted in Fig. 3, the step size being nor- 
malized by the standard deviation a. One observes that 
the distribution is slightly larger than that obtained if 
the x and y displacements were distributed according to 
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a gaussian centered on zero (solid line in Fig. 3). This 
observation is again very similar to the step size statistics 
observed in colloidal glasses [l(| . 

Another analogy with colloidal glasses is observed 
when studying the correlation between two successive cy- 
cles : the random motion seems to be anti correlated. 
When a particle moves in one direction, it has a great 
probability to move in the opposite direction the next 
step, which is reminiscent of a cage effect. This is shown 
in Fig. 3b where the component x' of the n + 1 th dis- 
placement along the direction of the n th displacement is 
plotted against the size Ar of the n th displacement. As 
in thermal colloidal glasses, x' is statistically negative, 
revealing anti correlation. 

The granular compaction system being not invariant 
in time, we study how the step size distribution varies 
during compaction. The different curves in Fig. 3 cor- 
respond to the distribution at different time in an ex- 
periment at 9 = 8.1° and show that the shape of the 
normalized distribution does not vary much with time. 
However, the mean value < Ar > varies during the com- 
paction as shown in Fig. 4. The mean step size < Ar > 
is plotted as a function of the volume fraction <3> during 
experiments carried out for three different shear angles 
6. In the three runs, the mean step displacement first 
decreases down to = 0.63 which corresponds roughly 
to 500 cycles and eventually reaches a constant value for 
dense enough packing. The asymptotic value depends 
on the shear amplitude and is smaller for small 0. An 
interesting observation is that once the packing is dense 
enough (<£ > 0.64), the mean step displacement is only a 
function of 6 independent of the history of the packing. 
This means that the same < Ar > is measured if the 
compaction is carried out at a constant or if the pack- 
ing is first compacted under a shear amplitude 9\ and 
then sheared at 6. This is shown in Fig. 5 where we plot 
< Ar > as a function of for many different experiments 
carried out at constant or non constant shear amplitude. 
The striking feature is the linear relation < Ar >= (30 
with (3 = Sjam/deg. The higher the shear amplitude, the 
larger the random displacement between two cycles. 

This linear relation between microscopic displacement 
and shear amplitude is reminiscent of the macroscopic 
observations of our previous study 7]. When a sudden 
change AO was imposed to the packing, a rapid (in less 
than 10 cycles) variation of volume fraction A<£ was mea- 
sured. A<3> was found to vary linearly with AO inde- 
pendently of the packing history for compacted enough 
packing. One can wonder if a link exists between the mi- 
croscopic random motion and the macroscopic behavior. 

In order to investigate this link, experiments are per- 
formed where the shear amplitude is discontinuously var- 
ied. A typical run is presented in Fig. 6. The packing 
is first compacted at a fixed amplitude = 8.4° for 10 4 
cycles and then is varied from 10.2° to 0.6° by step of 
1.2°, the angle being kept constant during 1000 cycles. 
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FIG. 4: Evolution of the mean step size < Ar > with <3> 
during compaction for = 2.7° (o); = 5.4° (x); = 8.4°(D). 
The average <> means over 40 identical slots in <E> and over 
particles. Fluctuations from one particle to another due to 
heterogeneities in the packing and force network are given by 
the error bar. 
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FIG. 5: Mean step size < Ar > as a function of shear angle 
0. The different symbols correspond to different experiments 
with kept constant or changed stepwise. 



The corresponding volume fraction variation is plotted 
in Fig 6a. As expected from our previous study 7], suc- 
cessive increasing steps in volume fraction are observed. 
The typical microscopic behavior of a particle during this 
experiment is presented in Fig 6b. The volume explored 
by the particle during its random motion successively 
shrinks when the shear amplitude decreases because the 
mean step size < Ar > decreases with (Fig. 5). How- 
ever, each time the shear angle changes, the other parti- 
cles below the test particle experience the same decrease 
in their exploration volume. The result is a net downward 
motion observed when the angle changes. The observed 
volume fraction variation thus seems to result from the 
change in the volume randomly explored by the particles. 
This becomes clear when looking at the displacement of 
all the particles during a change in (Fig. 6c). Thanks 
to the laser light visualization, we are able to follow each 
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FIG. 6: a) Volume fraction as a function of cycles when 
varies stepwise (see text), b) Corresponding trajectory of one 
particle. Changes in color correspond to changes in 0. c) 
Displacement field measured in the cell when changes from 
10.4° to 1.4°. 



ics of the system. 

In conclusion a simple scenario can be proposed for the 
compaction process. We have seen that the particles mo- 
tion can be divided in two kinds : random motion within 
confined volumes or cages, whose extent is directly pro- 
portional to the shear amplitude, and occasionally longer 
exploration corresponding to a change of cage. We think 
that the slow dynamics of compaction observed in exper- 
iment at constant amplitude may be attributed to the 
change in cages. This changes are irreversible and push 
the system towards more and more compact configura- 
tions. How cage deformations occur, what is the influ- 
ence of one cage rearrangement on the next ones are open 
questions. The cage changes certainly involve complex 
cooperative processes associated to modifications in the 
contact and force network which have not been studied 
here. However, we have shown that the rapid change of 
volume fraction observed when changing the amplitude 
is of different nature: it is simply related to the change of 
the cage size, without important structural changes. This 
means that this variation of volume is reversible and can 
be recovered by coming back to the previous amplitude 
of excitation. The existence of these two processes which 
affect differently the packing volume fraction, explains 
that memory effects can be observed. 

We thank F. Ratouchniak for his technical assistance. 



particle in a plane when the shear angle changes from 
10.4° to 1.4°. An homogeneous deformation is observed 
(Fig. 6c), as if when the angle decreases, all the particles 
deflate yielding a change in volume fraction. An estima- 
tion of the change in volume fraction can be made. As- 
suming that the random fluctuating motion of a particle 
occupies a volume of the order (d+ < Ar >) 3 , the volume 
fraction is <E> = ^ott^ — r — ttttt where <£n is a reference 

w (d+<Ar>(£/)) 

volume fraction. Knowing from Fig. 5 that < Ar >= (30, 
a change AO in the shear amplitude induces a volume 
fraction variation A<l> = — 3&o(3A0. A linear variation is 
predicted as observed in |g and the predicted coefficient 
(taking $ = 0.7) is -3<$> (3 = -0.56 10~ 3 deg- 1 hav- 
ing the right order of magnitude compared to the value 
— 1.5 10 -3 deg -1 measured in [3. 

All the results presented up to now are obtained for 
monodispersed glass beads. In this case, the compaction 
occurs towards crystallization |g . We have checked that 
the microscopic behavior does not qualitatively change 
when using a bidispersed material (mixture of 2 mm and 
3 mm beads) which do not crystallize. Cages are still 
observed and the sudden change of volume fraction when 
changing the shear amplitude |S| is correlated to a change 
of the mean random step size. The crystal nature of the 
packing then does not qualitatively influence the dynam- 
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